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Gauge-theory approach to describe Dirac fermioiis on a disclinated flexible membrane beyond 
the inextensional limit is formulated. The elastic membrane is considered as an embedding of 2D 
surface into R^. The disclination is incorporated through an SO{2) gauge vortex located at the 
("^ ' origin, which results in a metric with a conical singularity. A smoothing of the conical singularity 

^vj , is accounted for by replacing a disclinated rigid plane membrane with a hyperboloid of near-zero 

curvature pierced at the tip by the SO{2) vortex. The embedding parameters are chosen to match 
the solution to the von Karman equations. A homogeneous part of that solution is shown to stabilize 
_j , the theory. The modification of the Landau states and density of electronic states of the graphene 

■^Ui ■ membrane due to elasticity is discussed. 
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I. INTRODUCTION 



Is is now generally accepted that the low-lying electronic states in graphene can be accurately described by two- 
dimensional massless Dirac fermions on a plane ll|. In experiment, multiform CTaphene structures were observed 
C^ ' thus stimulating studies of Dirac fermions on curved graphene sheets (see, e.g., [2, |3|). This problem is markedly 
complicated when the curvature itself is generated by topological defects like disclinations. Indeed, a disclination is 
known in elasticity theory as a line defect which can be produced by "cut and glue" Volterra process, namely, by 
' O I inserting or removing a wedge of material with the following gluing of the dihedral sides. This immediately generates 
additional large elastic strains inside the crystal. For flexible membranes, however, there is a chance to screen out the 
strain field by buckling into a cone. The problem thus reduces to coupling Dirac spinors to a topologically non-trivial 
curved background. 

The topological lattice defects in graphene are pentagons that are equivalent to wedge disclinations. The first 

experimental observation of a pentagon at the apex of a cone was provided by An et al. in 4j . They used the scanning 

JT^ tunneling microscope to study the structure of a conical protuberance and found five bright spots at the apex of the 

^.f-\ nanocone. This was the first clear evidence that the pentagon is located at the apex. The bright spots indicate also 

^sj that there is an enhanced charge density localized at each carbon atom in the pentagon, which implies an increase 

C^ in the electronic density of states (DOS). This finding was confirmed by the numerical tight-binding calculations for 

ly-v , nanocones with different number of pentagons at the apex [5|| . Explicit manifestations of the topological effects in the 

f~^ ' electronic properties of disclinated rigid graphene surfaces have been discussed at length in the literature |6l-[l5| . 

f^ ' However, disclinations are the sources of long range elastic stresses which also modify the electronic states. The 

^~H strain-induced effects in the electronic structure of graphene are now a topic of intense research in strain engineering 

^ [ to produce basic elements for all-graphene electronics. In particular, designed strains can generate electron beam 

• ^ . collimation, confined states, quantum wires [la] as well as energy gaps and a zero-field quantum Hall effect [13] in 

^% ' graphene. To incorporate internal stresses due to topological defects one is supposed to go beyond the inextensional 

J^ , limit. Study of electronic properties of disclinated fiexible membranes beyond the inextensional limit is thus of 

- - I importance to address this new field of research. 

According to Volterra process disclination can be considered as a conical singularity like strings in cosmology. The 
relevant background is the curved spacetime where all the curvature is concentrated at the apex of the cone. The 
metric in the cylindrical coordinates is written as 

ds^ =dt^ +dz^ +dr^ +a^r^dip^. (1) 
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Here the parameter a is related to the angular sector that is removed or inserted to form the defect. In this case, any 
attempt to build a closed loop around the disclination line will result in a closure failure. The deficit angle is equal to 
27ra with a = \ — v where v is the Frank index, the basic topological characteristic of the disclination. The positive 
sign of V corresponds to the removing of a sector. In this case the space has positive curvature. Correspondingly, for 
negative v one has a cone of negative curvature. Eventually, the problem reduces to a Dirac equation in the curved 
spacetimc. By the change r ^ r" /a metric ([T]) can be brought into the form 

ds^ = dt^ + dz^ + r2("-i) (dr^ + r^dif^) (2) 

which describes a featureless cosmic string located at the origin. 

In spite of the elegant form of this approach, there is yet an important open question concerning the so-called 
core region of the defect. To the best of our knowledge, this problem was for the first time raised in cosmological 
models [I^, [l^ where long-range effects of cosmic string cores were studied. In geometric theory of defects, an 
influence of a disclination core on the localization of electrons and holes was investigated in [2(l|. In both cases, the 
tip of the conical singularity is replaced by a smooth cap while at large distances a typical cone with the deficit 
angle 27ra emerges. In cosmological models the curvature of an infinite string is confined within a cylinder of a small 
radius a (the core radius) that possesses a direct physical meaning: the string has a characteristic core radius given 
by a w 1/M, where M is a mass scale at which the string is formed. Accordingly, the relevant 2D piece of the metric 
can be taken in the form 

ds^ = P^{r/a)dr^ + aV^ V, (3) 

where the range of the angular coordinate is 99 € [0, 27r) and P{r / a) is a smooth monotonic function satisfying the 
conditions 

lim P{r/a) = i^, P{r/a) = 1, r > a. (4) 

For example, in [20| the so-called flower-pot model was considered when the curvature of the disclinated media is 
concentrated on a ring of radius a, which results in the formation of a "seam" on the cylinder. 

This approach is of interest in the description of linear defects with a certain interior structure (finite thickness of 
a string). However, the situation is more subtle for a disclination on an elastic 2D surface, e.g., on a graphene sheet. 
The specificity of elastic membranes lies in that they may change both their intrinsic and extrinsic geometries due 
to stretching and bending of the membranes, respectively. This occurs because a fiexible membrane can relieve the 
internal strain by buckling out of the plane into a cone. In the inextensional limit of the infinite rigidity, the stretching 
energy of a membrane of the radius R is proportional to the system size, i?^, whereas a buckled membrane has the 
energy proportional to log R. In this case a creation of a disclination by using the " cut and glue" process results in 
a true cone. This corresponds to a point-like disclination defect. In reality, however, a membrane possesses finite 
elasticity, so that one needs to go beyond the inextensional-limit approximation. 

As is known, the classical theory of elasticity introduces a characteristic velocity, but does not lead to a characteristic 
length. In defect theory the length scale is introduced phenomenologically through the core radius of the defect which 
appears as an external parameter. Accordingly, there is no room within the linear theory of elasticity for the description 
of the core region. In the rotationally symmetric case it represents a small disc with a certain radius tq. Physically, 
the parameter tq comes into play through the boundary conditions imposed on the stress tensor to insure that the 
internal strain inside a membrane be kept finite to prevent it from falling apart. It is this quantity tq that sets a 
relevant short-range length scale in the present problem (similar to a for cosmic strings) . 

In this paper, we attempt a variant of the self-consistent gauge-theory approach to account for both the smoothed 
apex and the topological characteristic of the disclination defect. To take into account finite elasticity, we go beyond 
the inextensional-limit approximation. In doing so, we invoke both the linear elasticity theory as well as a phenomeno- 
logical approach to account for the elastic deformations. Dynamical variables in our theory are the embeddings of 
a 2D elastic surface into R^ . Parameters of the embedding appear as the matter fields interacting with an external 
gauge potential that describes a disclination. Within the linear scheme the model recovers the von Karman equations 
for membranes with a disclination-induced source generated by the gauge field. We explicitly show that the elastic 
deformations emerge due to the defect. We compute the effective metric generated by the disclination. That metric 
is determined by the external gauge field and the dynamical embeddings. It contributes to a topologically nontrivial 
part of the spin connection that couples Dirac fcrmions to a 2D disclinated elastic manifold. 

The paper is organized as follows. Section II contains a brief account of a gauge-theory approach to describe 
fcrmions on a disclinated surface in the inextensional limit. To take into account elastic deformations, a dynamical 
theory beyond the inextensional-limit approximation is presented in Section III. Within that theory Dirac fermions 
couple to the underlying effective metric generated by the disclination defect in the presence of elasticity. In Section 
IV the derived approach is specified to study electronic properties of disclinated graphene. 



II. INFINITE RIGIDITY 

In the present Section we show that a dischnation on a 2D manifold can be thought of as arising due to an explicit 
breaking of the local rotational symmetry. This breaking can be enforced by a topologically nontrivial external gauge 
field that carries the dischnation "charge", v. To set the stage, we begin with a brief discussion of the case of an 
infinitely rigid membrane free to buckle into a cone in B? . This corresponds to the inextensional limit to be discussed 
from the elasticity theory point of view in the next Section. Here we just show that the conical-type singularity 
generated by the vortex results in a topological nontrivial contribution to the spin connection for Dirac fermions. Our 



discussion essentially follows that of our earlier short Communication [2l[ , except that it is extended now to include 
a dischnation dipole as well. 

A. Disclination vs gauge field 

Let a;" be a set of local coordinates on a Riemannian surface So- (Indices a,b,c, ... — 1,2 are tangent to So, whereas 
i,j, k, ... = 1, 2, 3 run over the basis of R^). To describe this we find it convenient to introduce an embedding Eg — >■ i?^ 
that can be realized in terms of a i?^-valued function K^^^Jx^ ,x'^). As the point (x^,a;^) is varied, vector _R(o) sweeps 
the surface Sg- This is nothing but a familiar two-parametric representation of surfaces in R^. In what follows the 
function R':qJx^,x'^) is chosen to specify an initial configuration Sq- Representation for the induced metrics follows 
immediately 

9ab = i9T.o)ab = daR(O) ' ^^^(o) . (5) 

Equation ([5]) is invariant under global SO{3) rotations of the vector R(^o)- To incorporate disclinations one should 
promote this invariance to a local one. To this end, consider the i??QN-bundle over Eq with the structure group SO{3). 

The so(3) valued one form {Wa ■ L)dx'^ serves as a connection one-form in the i??p>-bundle space over Eq, with Wa 

being the gauge potentials. Here L* £ so(3) are the generators of the group. By replacing in ^ ordinary derivatives 

daR(o) by the covariant ones Vq-R/q) — 9aR{o) + [Wa , R(o)], one arrives at the locally SO{3) invariant representation 
for the induced metric, 

9ab - .9i?(W>(")) - VaRiO) ■ Vbi?(0) = dailio) ' 4i?(0) + 9ai?(0) [V^f \ i?(o)] + 9f,^(o) [l^i°\ i?(0)] 

+ (W^i°^W>r)^(o) - (#i°)i?(o))(M>r^(o))- (6) 

Topological disclinations can then be considered as arising due to explicit breaking of the local rotational symmetry 
by a fixed topologically nontrivial gauge potential that generates a new metric as follows from Eq. ©. ^From now 
on a metric induced due to either the gauge field or elastic deformations will be denoted by gat to reserve the symbol 
9ab ^'-'^ ^^^ metric tensor on the undisturbed surface, Eq. 

In general case, a nonabelian gauge field W"*^") emerges to describe disclinated surface . However, throughout this 
paper we are primarily interested in the case Eq = R^ so that only the z-component of the gauge field matters. To 
illustrate this, consider a disclination defect placed at the origin of a plane that can be bent but cannot be stretched. 
We have Eq = i?^, so that (x^ = x, x^ = y) e R^. In this case W^°^'=^'^ = and WJ°^'=^ = W^°\ A singular 
vortex-like potential 



T^jo) = t/y/r^, W^°'> = -vxjr^, r = y^P^Ty^ 7^ (7) 

is supposed to describe a topological disclination with a strength v located at r = 0. This potential locally is a pure 
gauge. 

However, for any counter C encircling the origin one has 

W^'^Ux" = -2tti^ ^ 0. (8) 



/c 

Since the counter integral in Eq.® is a gauge invariant quantity, the field W\i cannot be gauged away to zero due 
to the topological obstruction. This is why that field is referred to as a topologically non-trivial one. A physically 



observable quantity associated with that gauge field is a nonzero flux, $ = — 27r;/, through an area bounded by the 
counter C. It does not depend on small continuous deformations of that area. This flux instead characterizes the gauge 
potential globally: it determines the first Chern characteristic class the gauge potential M^(°) belongs to. An electron 
encircling the origin naturally acquires a topological phase associated with that nontrivial flux; the Aharonov-Bohm 
phase which distinguishes the gauge potential T4^(°) from a trivial one. 

In the polar coordinates (r, ip) £ Rl^ a plane can be regarded as an embedding 

{r,(p) ^ {r cos (fi, rsm(p,0), < ?' < cx), < (/? < 27r. 

The gauge potential takes the form 

W,(°) = 0, VK^°) = -ly. (9) 

The components of the induced metric ([6]) can be easily read off 

Orr = 1, ffw = "^''^' 9rv = 9vr = 0, (10) 

Evidently, this is a metric of a cone (cf. (HJ) which at i^ = goes over to a flat one. 
In general case the potential 

^W;(-i...-«)(f) ^ _y-^Wf2^iL_|)!^ (11) 

A-^ r - n ^ 

where a,b — x, y, is supposed to describe N disclination defects with the strengths z^^*-' located at the points with the 
coordinates fi, i = 1, .., N. To specify to an important in applications case of the disclination dipole located on the 
X-axis, we put fi = {—L, 0), r2 — [L, 0) and i^*^^' = — 1/(^) = v. The above equation then yields 

where r\ = (a;±L)^ + J/^. Substituting this expression in Eq.® one can in principle compute the induced metric. To 
illustrate this, we work out an explicit representation of the metric induced by a disclination dipole in the asymptotic 
region Ljr <C 1. Far away from the dipole Eq. ([T2l) reduces in the polar coordinates to 

W'io) = -2j/ecos(^ + 0(e2), wW = — sin(^ + C'(£2)/r, £ := L/r < 1. (13) 

^ r 

Our representation of the induced metric ([5]) then gives 

ffrr = 1 + 0(2^), g<^y = r^(l-4iyecos(p + e'(£^)), .g^^) = .9^^ = 2i/ersin<p + re'(e^), (14) 

This metric describes the asymptotic representation of the line element for a screw dislocation (see, e.g., [22]), with 
the Burgers vector W — —AttvL being perpendicular to the defect region. The gauge-theory approach is thus seen 
to recover a well-known result: at large distances from the defect a disclination dipole can be thought of as a screw 
dislocation. 

B. Incorporating fermions 

Let us now move on to a problem of coupling fermions to a given disclination. As is known, the topologi cally 
nontrivial gauge field reasserts itself in the Dirac equation as a topologically nontrivial piece of the spin connection [23| . 
That part of the connection carries a topologically nontrivial flux that does not depend on smooth continuous changes 
of the underlying metric due to small elastic deformations. To incorporate fermions on the 2D background (Sq = 

R^^f ^ 0; gab {W^'^^)) we need a set of orthonormal frames {ea{W'^'^^)} which yield the same metric, 5q(, (W^*-"-*), 
related to each other by the local 5*0(2) rotation, 

e„->e^=A^e/3, Af e 50(2). 

It then follows that gah = e°e^Sap where e% is the zweibein, with the orthonormal frame indices being a, f3 = {1, 2}, 
and coordinate indices a, 6 = {1,2} (from now on we drop an explicit VF-dependence of the metric). As usual, to 



ensure that physically observed values be independent of a particular choice of the zweinbein fields, a local so(2)- 
valued gauge field w^ is to be introduced. The gauge field of the local 50(2) group is referred to as a spin connection. 
For the theory to be self-consistent, zweinbein fields must be chosen to be covariantly constant [24| : 

daet - Tl,eZ + K)^ef = 0, 
which determines the spin connection coefficients explicitly 

{uJar^^etDae^\ D, = da + Ta, (15) 

with Fa being the Fevi-Civita connection. The Dirac equation on a surface (Eo, fflb (^)) i^ written as 

il''el{da + na)i^ = E,l;, (16) 

with 

^a = i<^[7«,7/5] (17) 

being the spin connection in the spinor representation. In two space dimensions, the Dirac matrices can be chosen to 
be the Pauli matrices, 71 = — o"2, 72 = — cri. In the case under consideration Eq. ([15]) gives 

Hence, topologically nontrivial gauge field © results in a conical singularity of the spin connection. The flux 

ujfdif = 2ttv ^Q 

represents a "net" effect produced by a disclination on the moving electrons. We thus show that the gauge- field 
approach in the inextensional limit exactly coincides with the standard " cut-and-glue" procedure. 

However, a cone with a point-like apex is mathematical abstraction since in a real situation the media has a finite 
stiffness, which would inevitably result in a certain smearing of a conical singularity. Therefore, a proper description 
of the disclination implies a smooth deformation of the metric and at the same time one has to preserve a conical 
behavior far away from the origin. Although such a surface can effectively be approximated by a hyperboloid, we 
show now that one cannot incorporate finite elasticity into the theory by simply replacing a cone by a smooth surface 
that asymptotically approaches a cone far away from the origin. This would simply eliminate the defect. 

To illustrate this, consider an upper half of a hyperboloid as an embedding 

(^, 1^9) — >■ (asinh^cos(y9, asinh^siniy9,c cosh^), < ^ < cx),0 < (^ < 27r. (19) 

The components of the induced metric can be written as 

g^^ = a cosh £, + c sinh £_, g^^ = a sinh ^, g^^ = g^^ = 0, (20) 

which in view of (|15p gives for the spin connection coefficients 

a coshf 

The spin connection in the spinor 5*0(2) representation becomes 



, ,12 _ 21 _ r, 12 _ 21 



1- 



a;(C). (21) 



57;p — iujuz- (22) 

Since cl'(^) goes to zero as ^ ^^ a circulation of that field over a loop encircling the origin gives a flux which tends 
to zero as the counter shrinks to zero, 

lim ® ^^]ndii) = 0, 

where Og stands for a closed counter which encloses a small area ~ e^ around the origin. This equation implies that 
there is no a topologically nontrivial part in the flux. It is therefore clear that one should work out some other way to 
explicitly accommodate elastic deformations in the "cut-and-glue" procedure that would preserve a conical singularity 
at the origin. 



III. FINITE RIGIDITY 

In the present Section we formulate the hnear elasticity theory in terms of the embeddings, which proves convenient 
to incorporate defects. We thus arrive at the von Karman equations to describe elastic disclinated media away from 
the limit of infinite rigidity. Our consideration allows then to introduce an important notion of the defect core radius, 
tq. We assume that the linear theory works well in the region r > 7'o whereas it breaks down within the core, at 
r < Tq. We then briefly discuss the phenomenological approach discussed earlier to effectively account for the elastic 
properties of disclinated membranes |21| . Matching the linear von Karman theory at the boundary point r = rQ with 
the phenomenological approach enables us to represent the phenomenological parameters in terms of the elasticity 
and bending constants. 

A. Elastic surface 

Let us start by discussing the elastic properties of a 2£) manifold in the absence of defects. Under elastic deformations 
a surface Eq evolves into some other Riemannian surface E, which can be thought of as a diffeomorphic map, : Eq — > 
E. Again, we find it convenient to introduce the embedding E — >■ i?^ that can be realized in terms of a i?^-valued 
function i?*(x^,a;^), the point being that [25[ 

Rix):^rRio)^RioMx)], (23) 

where (/>* is a pullback of : Eq ^- E. The induced metrics becomes 

9at = (0 9^)at - (5s)cd^ ■d^-d^c-d^,Q^-d^-9aR- d,R, (24) 

where the set {0"} stands for local coordinates on E. The strain tensor is then determined to be 

^ab — gab ~ 9 ab ' 

The properties of a fiuctuating elastic surface are encoded in the action 

F^Fei+Ffi, (25) 

where F^i describes the elastic properties of the media, whereas Ffi stands for the Helfrich-Canham action to describe 
the energy of a free fluctuating surface. Explicitly the stretching energy is taken to be quadratic in the strain. 






(26) 



where trE = g^^.Eab: 5^"-' = '^s^llslb II ^^'^ summation over repeated indices is assumed. Here A and ^ are the 2D 
Lame coefficients. The stress tensor is then introduced to be 

aab = 2fiEab + XtrESab- (27) 

The Helfrich-Canham bending energy of a membrane depends on its mean curvature H and Gaussian curvature K 

MM. 



Fji^^ f y/^dx^dx^ H^ + ^ f Vg^dx^dx^ K 
2 Jeo 2 Jso 



(28) 



where k is a bare bending rigidity and kg is a Gaussian rigidity. H — g?Q-.Kab is the mean (extrinsic) curvature, and 
K — det g?I^^Kiic is referred to as the Gaussian (intrinsic) curvature. Here 

Kab^N- DaDbR (29) 

is the curvature tensor, and N is the unit normal to the surface 

\[d,R,d2R]\' 
The covariant derivative 

Da := da + Ta 

includes the Levi-Civita connection Fq. 



B. Disclinations in flexible membranes 

To incorporate disclinations originally distributed on Eq one needs to make in Ea.(|24p substitution 

daRiO) -^ Vai?(0) = daRiO) + [WJ,°\ R(^o)]- 

As a result, metric (/^j^ goes over to gab given by Eq.®. To derive equations of motion that follow from the Hamilton's 
principle of least action, 6F = 0, we first need to specify the embedding R^{x^,x^). In plane elasticity theory, an 
elastic deformation is represented by a displacement vector Ux{x,y),Uy{x,y). If a membrane is allowed to buckle out 
of the plane, we must add an extra function f{x,y) to describe the "deflection". We thus choose the embedding in 
the following way 

R{x\x^)^R^o) + U, (30) 

where i?(o) — (a;,y, 0) and U = {ux,Uy, f{x,y)) is a displacement of the (x,?/, 0) point under deformation. In the 
linear approximation we may omit the terms quadratic in displacements Ua as well as in the source strength v. The 
strain tensor then becomes 

Eab = daUb + dbUa + OafObf - £„„ W^i°^i?fo) " ^afeW^i^'^fo) + 0{u\ud!, W^) . (31) 

Introducing the Airy stress function x, 

fab = £ac^bddcddXj 

one eventually gets the following equations of motion 

KA'f = (d^yxMf) + {dlx){dlf) - 2{d,dyx){dxdyf), 

K-'A\ = {dxdyff - {dlf)(dlf) - eabdaWi''\ (32) 

where Ko — 4/i(A + ^)/(A + 2^) is the 2D Young's modulus. 

A single disclination located at the origin of a plane is described by the potential ([7]). This results in 

-eabdaWl°^ = ^/Alogr = 2^v5{r), 

so that Eqs.(I321) are exactly the von Karman equations for a defect in a flexible membrane [2^. Note, however, that 
the source term does not appear in ([32]) ad hoc but is rather generated by the gauge field due to a disclination. In 
case N disclinations with strengthes i^i, i = 1,2, ...,N are located at points r, one should use the gauge potential (|lip . 
The source term in Eg. ([5^ is then computed to be 



-CabdaWl-'—-"'' = 27^J2'^^Sir- n). 



Dynamically induced metric on E takes the form 

gab = 3ab + daUb + dbUa + dafdbf. (33) 

Both functions x and / or, equivalently, Ua and / contribute to it. Besides, function f{x, y) determines a shape of the 
emergent surface. To see this, let us reexamine the case of a disclination on a plane that can be bent but cannot be 
stretched considered in Section II (see Eg . (fTU|) ) . Within our approach it corresponds to the case of Kq — >■ cx). Since 
there is no in- plane stretching one may put Ua — Qva Egs. p2)) . The second von Karman equation takes the form 

[d.dyff - {dlf){dlf) = -2^v5{r), (34) 

whereas the induced metric becomes, 

gab = <5q6 + dafdbf. (35) 



Equation (|M|) possesses an obvious solution / = ±y/2h'r, which is a defining equation of a cone. The metric ((35 
coincides with that given by Eg. (1101) up to 0{v) order as it should in the linear approximation. 



C. Solution to the von Karman equations 



It is instructive to reveal a geometrical structure of Eq. (p4| . The embedding (pO| tells us that for an infinite stiffness 
the structure of the surface E is entirely determined by function z — f{x, y). In particular, the Gaussian curvature of 
the surface equals 

It is clear that a general solution to Ea.p4p must scale as / ~ y^. In the linear approximation one can therefore 
drop the /-dependent terms in the denominator in p6p. whereupon that equation takes the form 



K^2Triy6{r}. (37) 

This means that all the curvature is located in this case at the apex. Thus, to get the curvature spread over some 
finite area one needs to take into consideration elastic properties of a media. In other words, one needs to consider 
the whole set of the von Karman equations ((32|) at finite elasticity and bending constants. 

To this end, let us first rewrite the von Karman equations in the dimensionless form. Under the substitutions 
X -^ XK, r — > fro and / — )■ fro where rg is yet unspecified parameter with the dimension of length, those equations 
become 

AV = id^yX)id^f) + idh)id^yf)-nd.dyxKd,dyf), 

eA^X = 27rv6{f)-K, (38) 

where the parameter e — K/{KQr'^) and the Gaussian curvature K — det{dadbf)- All the functions as well as the 

1 /^ 
coordinates entering these equations are now dimensionless. Since in 2D one gets [Kq ] = [i^] — E, the parameter 

e is dimensionless as well. The inextensional limit amounts to that of e — >■ 0. Let us denote the solutions to the von 

Karman equations in this limit by /o and Xo- The second line in Eas. (|38l) becomes 



{dlfoKd'yfo) ~ {d.dyhf = 27ny6{r} 



which has an obvious solution /o = ±v2i^r that describes as already mentioned a true cone. Inserting /o into the 
first equation of (|55|) gives 

1 = (x^C + y'd'yy + 2xydly)xo. (39) 

This equation possesses an obvious inhomogeneous solution xo — — logr discussed at length in [28| . However, there 
is also a nontrivial solution of the corresponding homogeneous equation missed in [28|. A general solution to Eq. pop 
must include a homogeneous term and reads 

Xo = -logr + (jr. (40) 



Here q is an arbitrary constant. This becomes evident upon rewriting Eq. (|39p in the polar coordinates, 1 = r'^d^.^xo- 
The homogeneous part of the solution ()40[) turns out to be of the utmost importance in stabilizing the theory and 
bringing out the physical meaning of the parameter rg . 

To see this, let us for a moment restore an explicit vq dependence in Eq. (|40p : 



Xo = -log— + g— . (41) 

Consider an elastic thin disk of the radius R with a disclination sitting at the origin. Because of rotational symmetry, 

anp = -dr{-d^x) 
vanishes identically for (|4T|). The radial component of the strain tensor 

cr„. = -drX + —dl^X 



yields, however 



1 fl , . 



It is clear that (Trrir = R) vanishes as i? — >■ (X). It is singular, however, at the origin. Moreover, if we ignored the 
homogeneous solution in (j42p we would run into a serious problem. Namely, a typical way to remove the singularity is 
to delete a small disk of material around the origin. It is easy to see that for q — the strain tensor arr behaves like 
— 1/r^ at the boundary of the excised disk. This is clearly physically unacceptable, which signals the instability of the 
theory (see discussion in [23 )• It should be stressed that the general solution (|4T|) allows us to avoid this difficulty. 
To show this, let us delete a small disk of the radius tq around the origin and require that cTrrir = ro) = ctq. This 
yields q — 1 + (JgrQ. If one requires vanishing stresses at the inner boundary, one should put ctq = 0. Therefore, the 
parameter tq characterizes the core region of the disclination. 

It is also important to note that the homogeneous term in the solution ([1^ significantly affects the stretching 
energy of the membrane, 

Es = ^f d'r{y\r. (43) 



2Kq 



r>ra 



Since V^ logr = 2Trd{r) the stretching energy (|43p for the stress function P^ with q — becomes an identical zero. 
This is exactly the conclusion reached in [2g]. However, this is not physically appropriate, since this result should 
follow only in the limit Kq -> cx). If we instead compute (|43|) at (/ 7^ 0, we will get 

Es oc -—3 log — = enq^ log — , 44 

which indeed vanishes as K^ — >■ cx). It should be stressed that the same logarithmic behavior has the bending energy, 
which is written as [28i] 

£;b = 27ri/Klog— . (45) 

ro 

At g ^ 1 one has Eg ^ eEi, as it should be at small e. Because of the fact that the entropy in the Kosterlitz-Thouless 
argument also increases logarithmically, this result provides an interesting possibility of disclination-mediated phase 
transitions that might be realized in 2D elastic membranes (see, e.g., |29|). 

D. Phenomenological approach vs von Karman equations 

Let us now turn back to the dimensionless set of the von Karman equations ([55)1 . We seek a general solution in the 
form 

00 00 

n=0 n=0 

where we have already found /o and xo- These series are supposed to converge, provided the linear approximation 
is valid. Inserting pS)) back into ([55]) results in a set of the self-consistent coupled equations to determine step by 
step the functions /„ and x„. Technically, those equations for n > 1 turn out to be quite complicated. Their analysis 
will be given elsewhere. To get some insight, we invoke instead a sort of phenomenological approach to effectively 
incorporate elasticity in the Dirac equation discussed earlier in [2l[. We match the information that follows from 
our approach with that provided by the the von Karman equations. This enables us to explicitly determine the 
dependence of the phenomenological parameters on the elasticity and bending constants to analyze the electronic 
properties of elastic graphene. 

To begin with, let us briefly recall the phenomenological theory [2l|. As was already mentioned, a rigid plane 
pierced by a vortex results in a conical singularity. Let us now assume that the membrane possesses a small finite 
elasticity. In that case the vortex will produce the singularity at the origin as before and, additionally, it causes 
the medium to respond by smoothing the conical shape due to elasticity. We suggest that both effects can be taken 
into account by placing the vortex on the tip of a hyperboloid of a near-zero curvature. That sort of hyperboloid 
is supposed to effectively emerge as a response of the elastic plane to a disturbance caused by the defect at large 
distances, r > tq. The parameters of the hyperboloid must fulfil some natural requirements to be formulated shortly. 
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In this way we arrive at the effective metric that takes into account a response of the elastic media to the disturbance 
caused by the defect. 

Explicitly, we employ the embedding (|19p . with dimensionless (scaled by tq) parameters a and c that appear as the 
phenomenological parameters of the theory. The gauge field © represents the vortex at the origin. It induces the 
following metric: 

g^^ = a cosh ^ + c sinh ^, g^^ = a a sinh ^, g^^ = g^^p = 0, (47) 

where a = 1 — i/ and v is the declination charge. At t/ = this metric reduces to that of a true hyperboloid given 
by (1^0]) . Comparing the embedding IHM with Eq. (150)) tells us that c must scale as ^/ly, which seems natural: buckling 
is induced by the source. Let us further suggest that the elasticity coefficients are accumulated in parameter a, the 
inextensional limit corresponding to a — >■ cx). In other words, we should get that a — > oo as Kq — > oo. Let us check 
these two assumptions against the metric (|47p . We see that at i^ = it becomes 

344 = a^ cosh^ ^, g^^ = a^ sinh^ £,, g^^ = g^^ = 0. (48) 

Upon introducing a new variable r = a sinh ^ one gets 

df 

9rr = gii{r){^Y = I, giflip^r'^, g^pr^ gr<p^O, (49) 

which is nothing but the metric of a plane. This result is quite reasonable, since putting v — should result in the 
trivial solution to the von Karman equations, / = X — 0- In the second case of interest, we have v ^ and a 3> 1. 
With the help of the above mentioned change of variables, we get 

C V 
5rr = 1 + (-)^^-- J7 5¥'¥' = "^''^' g^pr ^ g^p = G- (50) 

We see that for large enough a, which in view of our assumption implies large K^^ the elasticity properties come into 
play through the single dimensionless parameter 77 — (? jc? ^1. In particular, in the inextensional limit a — > cxd, 
Eq. ([50| reduces to the metric of a cone. 

We now argue that the desired form of 77 follows directly from the von Karman equations. To this end we suppose 
that there exist two dimensionless functions x ^-nd / = /o + / that fulfil Eqs. psp and at the same time account 
for the embedding (fT9| as well as metric (|47| . We have explicitly singled out the /o dependence to eliminate the 
(5- function source. We are interested in the second "dynamical" equation in Eq. ([55)) . It takes the form 

eA^X = -KifoJ), r>l. (51) 

To proceed we simply replace the Gaussian curvature in this equation with that of the hyperboloid 

^ " (a2 + r2(l+?7))2- (^^) 

We recall that coordinate r as well as parameters a and b are now dimensionless. Equating then the both sides of 
(I5T]) at the core boundary r = 1 yields 

T]^ = kc^e + 0{e^). (53) 

In getting this we have assumed that a ^ 1. Here k — A'^x \r=i- To restore an explicit ro-dependence one needs to 
make the substitution a — >• a/ro, c —J- c/ro. This yields 

7,2 = fc(c/ro)2e. (54) 

It seems reasonable to assume that c ^ \/iyro, which finally gives 

r] ^ ^Jvl. (55) 

Eq. (|55)) qualitatively relates the characteristics of elastic media to the parameters of the embedding. It must be 
stressed that the smooth hyperboloid parametrized by EQ. (fT9)) at t^ = is not a surface that replaces a rigid cone for 
finite elasticity. Only when that hyperboloid of the near-zero curvature is pierced by the gauge flux, the smoothed 
replacement of the cone emerges. By construction a ^ tq and c '-^ rg, which is consistent with 77 <C 1. This takes care 
of a small stretching in the system. 
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Turning back to fermions, we see that the metric (j47p generates the spin connection term 

aa cosh^ 



12 21 



1 



/.9« 



^a(0- (56) 



Since ajQ,(^) — > 1 — a as f — > 0, it in contrast with (j2ip contains a topologicaUy nontrivial part that gives rise to a 
fixed flux, 



lim (p uj^^dw = 2tiv. 



V 



We thus get the smoothed apex, the cone-like asymptotic at large distances and the unremovable conical singularity 
at the disclination line. It is known that in case a spin connection contains an 5*0(2) piece with nontrivial flux, 
that field cannot be eliminated under any smooth deformation of the underlying metric (see, e.g., [2J]). Within our 
approach this simply means that a nontrivial contribution to the spin connection which comes from the topological 
gauge field survives any smooth elastic deformations of the media. 

IV. FLEXIBLE GRAPHENE 

In this section, we apply the developed approach to describe the electronic properties of graphene with pentagonal 
defects. Indeed, elastic characteristics of graphene are well-fitted to our theory. The estimated bending rigidity of 
graphene lies in the range of 1-2 eV and anyway does not exceed the value of 2.5 eV (see, e.g., |3C)| - [33 |). At the same 
time, the lower range value of Kqt^ is approximately given by 20 eV at rg = ao with ao being the interatomic spacing 
in graphene lattice [30|. Therefore, the parameter e is estimated as e < 0.1 thus justifying an applicability of the 
elasticity-induced perturbation scheme to graphene. 

A. Uniform magnetic field: Landau states 

The Dirac equation on a surface S in the presence of the gauge field a^ and the external magnetic field with the 
vector potential A\, is written as 

*7"ej'[V6-za6-afc]V^ = £;V, (57) 

where Vj, = dh + ^h- The effective abelian gauge field a\y is responsible for valley mixing since K and K' points 
become inequivalent in the presence of the pentagonal defect (see, e.g., [l| for detail). The energy in (|57l) is measured 
from the Fermi level. 

On a surface of the hyperboloid the Dirac operator reads 

le-(;|^ + ,rd^(^a,-i..„(e)-HM j' ^ ^ 

where 

ri<p = a^ + A^. (59) 

It can be verified that D = & . The substitutions 






i = 0,±l,..., (60) 



and 



V" = V'Vsinh^, (61) 

reduce the Dirac equation _DV = ^V' to the form 



(j -h 1/2 - fly + A^) 



d^u - "■■''' ^ ^\jco\h^i + r)u^Ev, (62) 
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, a + 1/2- a^+A^) I ^ . ~ . 

— OfW ^^ ^— A/coth ^ + riv = Eu, (63) 

a * 

where i? = JIJ^E. Explicitly, a^p = ±{N/A + M/3) for even nuniber of pentagons and a^p = ±N/4 for odd number 
of pentagons. Here N is the number of pentagonal defects and M takes the values M = —1, 0, 1 depending upon the 
arrangement of pentagons (see [33] for detail). The uniform external magnetic field B is chosen to be pointed in the 
z-direction so that A = B {y, —x, 0) /2. In (^, ip) coordinates the components of Ai, are written as 

A<p ^-^sinh^C; Aj^O, (64) 

where $ = ba'^/2 and b = eB/hc. Let us start with the analysis of the electron state at the Fermi level (so-called 
zero-energy mode) . In this case, we put E=0 in (|62p and ([M]) . The exact solution is found to be 

~ ,^^ ^/A , , ^.fci+ni M-fcoshA"^ / -lAcoshA 
uoiO=CiA + kcoshO^+-^'' i———^\ expl ^J, (65) 

~fc\ n'(/K^u ^f^-fcJ-l? /^^±cosh^y /$Acoshe\ .„„. 

vo{£,)=C' {/\ + kcoshS,) ( sinh£ j '^'^^ 1 2 1' ^ -^ 

where 



fc= ^1 + ry; A = A(e) = yi + fc^sinh^^, ~j = {j + 1/2- a^)/a, $ = $/a, 

and C and C" are the normalization factors. Evidently, the only component mq becomes normalizable. In the 
inextentional limit (77 -^ 0) one obtains 

uo(r) oc rj'-i/2 exp (-^) , (67) 

where r — asinh^. This result agrees with that of [33] • 

Let us consider the Landau states. To this end, we develop the perturbation scheme using 77 ~ i/e as the perturbation 
parameter. Indeed, according to (I50p the hyperboloid can at small -q be considered as a local perturbation of the cone 
metric. Notice that a similar procedure was used in the description of electronic states in spheroidal fuUerenes [34l. 1351] 
where the spheroid was considered as a slightly elliptically deformed sphere. At the same time, the Landau states 
on the graphene cone were already studied in detail in [Sj]. Therefore, one can use the unperturbed solutions found 
there. 



The Dirac operator is written as 



V = Vo+vA, (68) 



where 

a cosh f \ 2 sinh £ / a sinh £ v / 



is the Dirac operator on the true cone. It is important to note here that generally the operator Vi is not Hermitian 
on a cone and must be extended to a Hermitian one (see, e.g., the discussion in [341 1351]). The result is 

23i = -^^i^(j + |.sinh^£). 
2acosh^£V-' ^J 

It is convenient to square the Dirac operator, 

t>^^t)l+ r;f , (69) 

where T — 'Do'Di + 'Di'Dq and the quadratic in rj term is omitted. Explicitly, 
1 



a^ cosh £ 



j + ^smh'^A^ + ^ (j ( % 1) +I>sinh2£ 



2 V \cosh £ / Vcosh £ 



(70) 
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For r] = both the unperturbed wave functions ipjn and the Landau energy levels were obtained in [3^ . Two families 
of solutions were found. We restrict our consideration here to the first family where j > 0. In this case, the energy of 
the so-called bulk levels reads [3J| i?° = ±-\/2n, n = 0, 1, 2..., where the energy is measured in units of hvp/ls with 
the magnetic length l^ ~ yfic/eB. 

Let us calculate the matrix element of the perturbation F — (^j„|F|^j„). Our analysis shows that the perturbation 
does not influence the zero energy level in the first order in 7]. Since the resulting expressions are rather involved, we 
calculate the first energy level numerically. The first Landau level is found to be slightly shifted due to the elastic 
contribution, 

^.Li-±V2±ry^~±V2±0.3r;, ,7 > 0, (71) 

where the ± sign corresponds to the conduction and valence band, respectively. Note that a similar shift of the first 
Landau level follows for the second {j < 0) family of the solutions referred to as the apical states [33] . 

In graphene there are many allowed transitions due to the presence of two electron bands, the conduction and the 
valence band, and the transitions have the energies 

Ai+i,„ = ^[V2^^^TT)-CV2^], (72) 

where ^ = ± denote the intraband and interband transitions, respectively [36|. As is seen from Eq. (j7ip the transition 
between the ground energy level and the first one in the presence of elasticity is modified to become 

Ai.o = ^[y2 + 0.3?7]. (73) 

IB 

The experimentally observed cyclotron resonancelike and electron-positron-like transitions are in a good agreement 
with the theoretical expectations of a single-particle model of Dirac fermions in graphene [37[ . They produce a very 
accurate value for vp, the velocity of electrons in graphene. In our case, the Fermi velocity becomes slightly renor- 
malized due to the elasticity effects induced by the phenomenological parameter rj : Avp/vp = 0.3?7/'\/2. Therefore, 
it would be interesting to perform analogous experiments with graphite cones instead of graphene to determine the 
modification of vp. 

B. Zero magnetic field: Electronic density of states 

Let us consider now the case of zero magnetic field, B = 0. This markedly changes the situation with the zero-mode 
states. Indeed, in the absence of the cut-off exponent in (|65|) and (|66p the normalization conditions are found to be 
— 1/2 < j < —l/2k for Mo(C) and l/2fc < j < 1/2 for uo(C)- As a result, at small ry which is of interest here, there are 
no normalized solutions. This means that smoothing has no marked effect on the existence of zero modes. As stated 
above, one of the modes (either m(^) or u(^)) becomes in the presence of the uniform magnetic field directed along 
the 2;-axis nonmailable and there exists a true zero mode. Therefore, one can expect a "switching-like" effect driven 
by the magnetic field. 

An interesting question is how a hyperboloid geometry infiuences the density of states near the Fermi energy in the 
vicinity of the pentagonal defects. We will be interested in the DOS in a small ring < r < 6 around the defects. 
First, we need to find corrections to the wave functions for the hyperboloid geometry. Following the perturbation 
scheme with 77 being a small parameter, one can write i/jn = V'n + V'n + ••■ with -0^ being the solutions for 77 = which 
are the ordinary Bessel functions (see [S|)- Here ip^ oc J^m^mn'tpm/ '^^ is a perturbative part with Ak = tt/R. The 
matrix element of the perturbation reads Fm„ — (V'?m|r|V'?n); where the perturbation term takes the form 



J 



a^ cosh^ ^ 
Finally, the total DOS in the 5 disk is found to be 



^ 2 Uosh^e 



(74) 



D{E,5,n)o^D0S'^ +'qDOS'" oiDOS'^{l + r^\ — \5^), 77 > 0, (75) 



where 

i?05"(x|^|2"+M2("+l), 
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coincides with the total DOS near the defects for graphitic cones found in \§\, and n — j ± 1/2 is the index of the 
Bessel function. The term 

i:'OS"'cx |£;p("+i),52("+2), (76) 

comes from the elastic perturbation. Thus, we obtain an increase of DOS in the vicinity of the defect due to the 
finite elasticity. This behavior can be explained by the elasticity-induced contribution to the underlying metric which 
markedly deviates from that of a true cone near the apex. The modification of the DOS is more pronounced for 
more fiexible membranes. Notice that the bending rigidity of square graphene is found to grow with a size [31] while 
the Young's modulus depends on the temperature i32j . These findings give a possibility of experimental studies of 
the influence of elasticity on the electronic characteristics of disclinated graphene. As a promising material, one can 
consider the graphene monolayers obtained via chemical reduction of graphene oxide which have excellent mechanical 
properties including high bending flexibility and tensile strength '38'|. 

V. CONCLUSION 

In conclusion, we have presented an analytical approach to describe Dirac fermions on a flexible disclinated surface 
beyond the inextensional limit. The elastic membrane is considered as an embedding of a 2D surface into R^. 
The disclination is incorporated through a topologically nontrivial S0{2) gauge field that generates a metric with a 
conical singularity. A smoothing of the conical singularity is accounted for by regarding the upper half of a disclinated 
two-sheet hyperboloid as an elasticity-induced embedding. Parameters of that embedding are chosen to match the 
solution to the von Karman equations. Away from but close to the inextensional limit, the Young's modulus Kq and 
bending rigidity k enter the theory through a dimensionless parameter, e = K/{Korl) ^ 1, where r^ sets the relevant 
short-range scale: a radius of the defect core. 

We argue that the homogeneous part of the solution to the von Karman equations for disclinated membrane is 
of the utmost importance in stabilizing the theory. This finding allows us to avoid the evident problem with the 
core radius mentioned in [23]. Surprisingly, the discovered homogeneous part of the solution has been missed so far. 
We calculate the stretching energy of the membrane to observe the logarithmic behavior with the membrane radius 
R similar to that of the bending energy. For rigid enough membranes we work out the self-consisting perturbation 
scheme with e being the small parameter. 

We apply a new approach to study the structure of the low energy electronic states of flexible graphene with a 
topological defect. We find a true zero-mode state in the presence of an external uniform magnetic field. The finite 
elasticity results in a smearing of the cone apex thus modifying the cone metric. The elasticity affects enter the 
observable quantities through the elasticity-induced phenomenological parameter 77 ~ ^/i^e. Qualitatively, the first 
Landau level in the conduction band is found to be shifted upwards, whereas the corresponding one in the valence 
band is shifted downwards. The total DOS near the tip is shown to increase due to elasticity. 

From the experimental point of view, it would be interesting to carry out experiments with elastic graphite cones to 
determine the modification of vp due to elastic deformations caused by a disclination defect. The elasticity effects are 
more significant for more flexible materials. Therefore graphene-based materials with excellent mechanical properties 
like those recently obtained via chemical reduction of graphene oxide in [38| would be of interest in these studies. 
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